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Abstract 

The micromechanical generalized method of cells model is employed for the pre- 
diction of the effective moduli of electro-magneto-thermo-elastic composites. These 
include the effective elastic, piezoelectric, piezomagnetic, dielectric, magnetic perme- 
ability, electromagnetic coupling moduli, as well as the effective thermal expansion 
coefficients and the associated pyroelectric and pyromagnetic constants. Results are 
given for fibrous and periodically bilaminated composites. 


1. Introduction 

The method of cells and its generalization, referred to as the generalized method of cells 
(GMC), is an approximate analytical micromechanical model which is capable of predicting 
the overall behavior of continuous and discontinuous (short-fiber and particulate) multiphase 
composites from the knowledge of the properties of the individual phases and their volume 
fractions. A review of the method has been recently updated in Aboudi (1996), where critical 
assessments of the method and its application by various researchers were outlined. As 
documented, many types of composites (e.g. thermoelastic, viscoelastic, nonlinear elastic and 
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viscoplastic) have been analyzed by the method, and the reliability of the predictions were 
demonstrated under many circumstances. For a detailed explanation of the two-dimensional 
case with four subcells see the recent textbook by Herakovich (1997). 

To provide a framework for GMC, a code referred to as MAC/GMC has been developed 
at NASA Glenn Research Center. Through MAC/GMC, various thermal, mechanical and 
thermomechanical load histories can be imposed, different integration algorithms can be 
selected, many different fiber architectures can be utilized and a variety of fiber and matrix 
constitutive models are available (see Arnold et al. (1999), a user guide for version 3.0). 

In an attempt to incorporate intelligent materials into GMC, this method has been ex- 
tended in a recent paper (Aboudi, 1998) and applied for the prediction of the effective elastic, 
piezoelectric, dielectric, pyroelectric and thermal expansion coefficients of multiphase com- 
posites with embedded piezoelectric phases. The predictive capability of GMC was assessed 
in this case by comparison with results of Dunn and Taya (1993) and Dunn (1993), and very 
good agreements were obtained. 

Piezoelectric materials exhibit interactions among mechanical, electric and thermal ef- 
fects. A detailed account on piezoelectric materials, actuators and their applications has 
been given in a recent monograph by Uchino (1997). A generalization of the analysis of 
electro-mechanical interaction effects would be the incorporation of magnetic field. Such a 
generalization would enable the analysis of composite systems that involve the interaction 
of electro-magneto-thermo-mechanical effects. The analytical modeling of such composites 
provides the opportunity to study the effect of controlling and altering the response of com- 
posite structures that consist of composite materials the phases of which are, in general, 
electro-magneto-thermo-elastic. Piezoelectric and piezomagnetic composites would be ob- 
tained merely as special cases of the general theory. A first step in developing such a hybrid 
analytical approach has been presented by Carman et al (1995) where a concentric cylinder 
model has been employed. More recently, Li and Dunn (1998) presented a micromechanical 
methodology for the analysis of electro-magneto-thermo-elastic composites, and the Mori- 
Tanaka mean field approach has been applied to generate the effective moduli of composites 
that consist of piezoelectric and piezomagnetic phases. 

In the present paper, the fully coupled electro-magneto-thermo-elastic constitutive equa- 
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tions are incorporated into GMC in order to establish the various types of effective moduli 
of multiphase intelligent composites whose constituents, in general, behave according to 
these equations. Piezoelectric and piezomagnetic phases can be obtained as special cases 
from these constitutive laws. Several comparisons of the predicted results with those pro- 
vided by Li and Dunn (1998) for two-phase composites are given. 

2. Basic Equations 

The constitutive equations that govern the interaction of elastic, electric, magnetic and 
thermal fields in a electro-magneto-thermo-elastic medium relate the stresses cr^, strains 
6ij, electric field E t , magnetic field H z and temperature deviation AT (from a reference 
temperature T 0 ) as follows (see Li and Dunn (1998), for example): 

&ij Eijisi^ki qijkHk A-ij AT j, h,l = 1, ...,4 (1) 

where Cijki, e^k, qijk and A y denote the fourth order elastic stiffness tensor, the third order 
piezoelectric tensor, the third order piezomagnetic tensor, and the second order thermal 
stress tensor of the material, respectively. 

In addition, the electric displacement vector, A, is expressed in terms of the strain, 
electric field, magnetic field and temperature in the form 


Di — eikitki + K %kEk + o-ikHk — PiAT (2) 

where K ik , a ik , and are the second order dielectric tensor, the second order magnetoelectric 
coefficient tensor, and the pyroelectric vector, respectively. 

Finally, the magnetic flux density vector, is given in terms of the mechanical, 
electric, magnetic fields and temperature by 


Bi — qiki^-ki + a tkE k + p ik H k — rriiAT (3) 

where p ik and rn t are the second order magnetic permeability tensor and the pyromagnetic 
vector, respectively. These tensors satisfy certain symmetry conditions that are given by 
Nye (1957). 
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Let the vectors X and Y be defined as follows: 


X— [eu, 622 , ^33>2e23, 2ei3, 2e 12 , — E\, — E 2 , — E 3 , — Hi, — H 2 , — H 3 ] (4) 

Y = [(T 11, <722, °33) ^23) ^13 j & 12) D 2 , D 3 , Bi, B 2 , B3] ( 5 ) 

Consequently, equations (l)-(3) can be written in the following compact matrix form: 


Y = ZX - TAT 


(6) 


where the square 12th order symmetric matrix of coefficients Z has the following form 


and 



( 7 ) 


A 


r = P 


(8) 


{ m J 

In equation (7), the square matrix C of the 6th-order represents the 4th order st iffn ess tensor, 
e* denotes the transpose of the rectangular 3 by 6 matrix e that represents the corresponding 
third order piezoelectric tensor, q* denotes the transpose of the rectangular 3 by 6 matrix 
q that represents the corresponding third order piezomagnetic tensor, k is a square matrix 
of order 3 that corresponds to the dielectric tensor, a is a square matrix of the 3rd order 
that represents the magnetoelectric coefficients, and fj. is a square matrix of the 3rd order 
that represents the magnetic permeability tensor. Finally, the 6th order vector A, and the 
two 3rd order vectors p and m in Eqn. (8) represent the thermal stresses, pyroelectric and 
pyromagnetic coefficients, respectively. 
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3. Micromechanical Analysis 


Consider a multiphase composite material in which some or all phases are modeled as electro- 
magneto-thermo-elastic materials. It is assumed that the composite possesses a periodic 
structure such that a repeating cell can be defined. In Fig. 1 , such a repeating cell is shown 
which consists of N a NpN-y rectangular parallelepiped subcells. The volume of each subcell 
is d a hpl y , where a, f3 , 7 are running indices a = 1, • • • , N a ; f3 = 1 , • • • , Np\ 7 = 1, • • • , N y in 
the three orthogonal directions, respectively. The volume of the repeating cell is dhl where 

N a Np 

d = £ da , h = Y, hp , l = L, (9) 

a=l p = 1 7 =1 

Any subcell can be filled in general by electro-magneto-thermo-elastic materials. Electro- 
magneto-thermo-elastic unidirectional long-fiber composites, short-fiber composites, porous 
materials, and laminated materials are obtained by a proper selection of the geometric di- 
mensions of the subcells and with an appropriate material fillings. 

The micromechanical model employs a first order expansion of the displacement in each 
subcell (a/3 7 ) in terms of the local coordinates located at the center of the 

subcell. 


u, ( q/37) = Wl {a ^ ] 


+ 4 Q) 4>i iafiy) 




+ 4 7) i>i {afh) 


1,2,3 (10) 


where are the displacement components at the center of the subcell, and 

and ^ (a/37) are the microvariables that characterize the linear dependence of the dis- 
placement on the local coordinates , x^ , x ^ } . In eqn.( 2 ) and the sequel, repeated 

Greek letters do not imply summation. 

The electric potential ^ a/3l> is also expanded linearly in terms of the local coordinates of 
the subcell: 

£ (a/?7) = £ 0 (q/?7) + x[ a) ^ (q/37) + xf ] £ 2 (q/J7) + 4 7) £ 3 (a/?7) (11) 

Similarly, the magnetic potential rf 0 ^ is expanded linearly in terms of the local coordi- 
nates of the subcell: 


V 


(a/?7) 


rj 0 + x[ } + x 2 r? 2 (a/37) + x 


i 7) 


( 12 ) 
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( 13 ) 


The components of the small strain tensor are given by 

i“ w = + a,<4“ w ) t,j = 1,2,3 

where d\ = d/dx[ a \ 82 = d/dx^ and O3 = d/dx^\ 

The components of the electrical field E in the subcell are obtained from the electric 
potential via 

= -d, &+* (14) 

The components of the magnetic field H\ apl) in the subcell are obtained from the mag- 
netic potential rj ' via 

H (°Py) = _ d . Jafh) (15) 

With the linear expansions of the displacements and the electric and magnetic potentials 
given by eqns.(10)-(12) the static equilibrium of the material within the subcell (a/3 7) is 
satisfied, namely 

do” = 0 (is) 

Furthermore, in the absence of volume charges: 


£>|“ m = 0 


(17) 


is satisfied, as well as 

Bj = 0 (18) 

The volume average of the stresses electric displacements D t and the magnetic flux 
density B l in the entire repeating cell (namely in the composite) is given by 


1 N a Np 

Y=57tE E E dahfl, 

afU a = 1 0=1 7=1 


(19) 


Similarly, the volume average of the strains e^, electric field components E % and magnetic 
field components Hi in the composite is given by 


1 N a N p 

* = -TT. E E E <*« h tr 

afU 0=1 0=1 7=1 


(20) 


In the framework of the GMC micromechanical model it is possible to establish a re- 
lationship between these two volume averages. For a piezoelectric composite such a rela- 
tionship has been obtained (Aboudi, 1998) by establishing the appropriate concentration 
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matrix that relates the electromechanical field in the subcell to the composite strain e* r 
The derivation of the concentration matrix is based on the application of the appropriate 
interfacial conditions between neighboring subcells as well as between repeating cells (to 
insure periodicity). In the present electro-magneto-thermo-elastic analysis these interfacial 
conditions are: (1) continuity of displacements, (2) continuity of tractions, (3) continuity of 
the electric potential, (4) continuity of normal electric displacements, (5) continuity of the 
magnetic potential, and, (6) continuity of the normal magnetic flux density. By generaliz- 
ing the isothermal analysis (with AT = 0) that was implemented in the piezoelectric case 
to the present electro-magneto-elastic composite, one can establishes the following relation 
between the local electro-magneto-elastic field in the subcell X^ a/5 ^ and the average external 
(composite) macro field X in the form 

X (Q/?7) = A X (21) 

where is the electro-magneto-elastic concentration matrix associated with subcell 

Mt)- 

Substitution of eqn.(21) into (6) (with AT = 0) yields 

Y( q / 37) _ 2( Q ^7)^(«h7) X (22) 


Consequently, in conjunction with the averaging procedure given by (19), the following 
effective isothermal constitutive relations of the electro-magneto-elastic composite can be 
established 

Y = Z* X (23) 


where the effective elastic stiffness, piezoelectric, piezomagnetic, dielectric, magnetic perme- 
ability and electromagnetic coefficients matrix Z* of the composite is given by 

Np N 

E 

0=1 7=1 

The structure of the square 12th order symmetric matrix Z* is of the form 


Y d a hp L, Z {a ^> A ( ^ 7) 


(24) 


1 


N a 


= — Y 

dM Q= i 


C* e* 4 q' 



(25) 
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where C*, e*, q*, k* , a* and (i* are the effective elastic stiffness, piezoelectric, piezomagnetic, 
dielectric, magnetic permeability and electromagnetic coefficients, respectively. 

In order to incorporate the thermal effects in the composite, we utilize Levin’s (1967) 
result to establish the effective thermal stress tensor, A? , pyroelectric coefficients, p*, and 
pyromagnetic coefficients, m*. This approach has been also followed by Dunn (1998) to 
establish the requested effective thermal moduli. To this end let us define the following 
vector of thermal stresses, pyroelectric and pyromagnetic coefficients material constants in 
the subcell (afiy): 

T (q/37) = [Ai, A 2 ,A 3 , A 4 , A 5 ,A 6 ,p 1 ,p2,P3,mi,m2,m 3 ] (a/37) (26) 

The corresponding global (effective) vector is defined by 


r* = [AJ , A£ , Ag A 4 , A£ , Ag , p* , pj , P 3 , mj , , m^] 


(27) 


According to Levin’s result, the relation between and T* is given in terms of the 

electro-magneto-mechanical concentration matrices A^ 1 ' 1 as follows. 

I N a N P ^7 

r* = — g g E d a hp A 1 ^ (28) 

ant Q= 1 7=1 

where A^ a/3 '^ is the transpose of A^“^ 7 \ The above relation provides the effective thermal 
stress A*, pyroelectric p* and pyromagnetic m* vectors of the composite. 

Consequently, the final anisothermal micromechanically established constitutive law of 
the electro-magneto-thermo-elastic multiphase composite is given by 


Y = Z* X - T* AT 


(29) 


Finally, the effective coefficients of thermal expansion a* and the associated pyroelectric 
constants P* and the pyromagnetic constants M* of the composite can be assembled into 
the vector: 


rr = K, <*;, a;, a;, 0$, ag, Pi pi Pi Ml m*, m;\ 


(30) 


This vector is given by: 


si* = z* _ 1 r* 


(31) 
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It should be mentioned that a similar relation holds for the local quantities, namely for 
the the coefficients of thermal expansion and the associated pyroelectric and pyromagnetic 
constants of the materials filling the subcells. In such a case Z* and T* should be replaced 
by and r^ a/37 \ respectively. 


4. Applications 

Following Li and Dunn (1998), let us consider a composite consisting of a CoFe20 4 piezo- 
magnetic matrix reinforced by BaTi0 3 piezoelectric material. Both phases are transversely 
isotropic with the axis of symmetry oriented in the 3-direction. The independent material 
constants of these constituents are given in Tables 1-3 (Li and Dunn, 1998). It should be 
noted that in both materials the electromagnetic coefficients are zero, i.e., a l3 =0. 


Table 1. Elastic material properties. 


Material 

C n (GPa ) 

C l2 (GPa) 

Cn(GPa) 

C 33 (GPa) 

C 44 (GP a) 

BaTi0 3 

166 

77 

78 

162 

43 

CoFe-zO^ 

286 

173 

170 

269.5 

45.3 


Table 2. Electric material properties. 


Material 

eis (C/m 2 ) 

e 31 (C/m 2 ) 

633 {C/m 2 ) 

kh(10- 9 C/Fto) K 33 (10~ 9 C/Vm) 

BaTi0 3 

11.6 

- 4.4 

18.6 

11.2 

12.6 

CoFe20 4 

0 

0 

0 

0.08 

0.093 


Table 3. Magnetic material properties. 


Material 

qis(N/Am) 

q 3 \{N/Am) 

q 33 (N/Am) 

± * 

fj, n (10~ 6 Ns 2 /C 2 ) 

^33(10 ~ 6 Ns 2 /C 2 ) 

BaTi0 3 

0 

0 

0 

5 

10 

CoFe20 4 

550 

580.3 

699.7 

-590 

157 


Consider first a fibrous composite in which the BaTiO 3 continuous fibers are oriented in 
the 3-direction. The volume fraction of the fibers is denoted by vj. In this case the GMC 
approach can be employed to generate the effective elastic, dielectric, magnetic permeability, 
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piezoelectric, piezomagnatic, and electromagnetic coupling moduli of the fibrous composite 
for 0 < Vf < 1. To this end, the composite is modeled by the simplest configuration in 
which two subcells exist in every one of the orthogonal directions, see Fig. 1. These moduli 
were also generated by Li and Dunn (1998) who employed the micromechanical model of 
Mori-Tanaka (1993), MT. Figures 2 through 7 show a comparison between these effective 
moduli predicted by GMC and MT. It is readily seen that in some cases the predictions of 
both methods coincide, while in other cases some deviations are observed. It should be noted 
that both GMC and MT are approximate models (for some shortcomings of MT model, see 
Ferrari (1991)), and it is not possible to determine which one is more accurate. 

Li and Dunn (1998) also generated the effective moduli of a bilaminated composite that 
consists of alternate layers of BaTiO 3 and C'oFe 2 0 4 , in which the layering is in the 3- 
direction. Figs. 8-13 exhibit the corresponding effective moduli for this type of composite 
as predicted by GMC and MT. In this case the predictions of the two methods coincide for 
all values of volume fraction Vf of the BaTiO 3 phase. 

As mentioned by Li and Dimn (1998), the most interesting behavior is the overall electro- 
magnetic coupling effect that is present in the composite, but not in either of the individual 
phases. This effect is shown in Figs. 8 and 13 for fibrous and bilaminated composites, respec- 
tively. As discussed by these authors, the different order of magnitudes of the coefficients 
ajj and O33 exhibited in both figures is attributed to the difference in magnetic constants 
between the piezoelectric BaTiOz and piezomagnetic CoFe 2 0 4 phases. Further discussion 
of the various effective moduli behavior has been given by Li and Dunn (1998). 

5. Conclusions 

The fully coupled electro-magneto-thermo-elastic constitutive law has been incorporated 
into the GMC approach, thus providing the ability to model advanced composites that are 
sensitive to electric and magnetic fields. Since the GMC is a multiphase micromechanical 
model, it can be easily employed to investigate a wide range of effects, as has been shown by 
several investigators. For example, one can investigate the effect of inelasticity of the host 
material, the effect of weak bonding between the phases, the effect of fiber distribution and 
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architecture on the overall response of composites, and since GMC is a constitutive law, it can 
be utilized by a structural analysis package to investigate the behavior structures composed 
of these corresponding composite materials. Now with the electro-magneto generalization to 
the GMC outlined, herein, an entirely new class of materials can be analyzed and designed 
with this powerful approach. 
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Figure 1 . — A repeating cell in GMC consisting of N ai 
A/p and N y subcells in the 1 , 2 and 3 directions, 
respectively. 


Figure 3. — Effective dielectric moduli of fibrous com- 
posite against the volume fraction of Ba770 3 . The 
predictions of GMC and MT coincide. 



Figure 2. — Effective elastic moduli of fibrous compos- 
ite against the volume fraction of Ba770 3 . The pre- 
dictions of GMC and MT coincide. 



Figure 4. — Comparison between GMC and MT pre- 
diction of the effective magnetic permeability 
moduli of fibrous composite against the volume 
fraction of Ba770 3 . 
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Figure 5. — Effective piezoelectric moduli of fibrous 
composite against the volume fraction of Ba770 3 . 
The predictions of GMS and MT coincide. 



Figure 6. — Comparison between GMC and MT predic- 
tion of the effective piezomagnetic moduli of fibrous 
composite against the volume fraction of BaTiO 3 . 



Figure 7. — Comparison between GMC and MT predic- 
tion of the effective electromagnetic coupling moduli 
of fibrous composite against the volume fraction of 
Ba770 3 . 



Figure 8. — Effective elastic moduli of bilaminated 
composite against the volume fraction of Ba7/0 3 . 
The predictions of GMC and MT coincide. 
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(KHNs^C 2 ) ^ (1 0 -9 C/Vm) 



igure 9. — Effective dielectric moduli of bilami- Figure 1 1 . — Effective piezoelectric moduli of bilami- 

nated composite against the volume fraction of nated composite against the volume fraction of 

Ba770 3 . The predictions of GMC and MT coincide. Ba770 3 . The predictions of GMC and MT coincide. 



Figure 1 0.— Effective magnetic permeability moduli Figure 1 2.— Effective piezomagnetic moduli of bilami- 

of bilaminated composite against the volume nated composite against the volume fraction of 

fraction of Ba7/0 3 . The predictions of GMC and Ba770 3 . The predictions of GMC and MT coincide. 

MT coincide. 
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Figure 1 3. — Effective electromagnetic coupling moduli of 
biiaminated composite against the volume fraction of 
BaTiO$. The predictions of GMC and MT coincide. 
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permeability, electromagnetic coupling moduli, as well as the effective thermal expansion coefficients and the associated 
pyroelectric and pyromagnetic constants. Results are given for fibrous and periodically bilaminated composites. 
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